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General Relativity obeys the three equivalence principles, the "weak" one (all test bodies fall the 
same way in a given gravitational held), the "Einstein" one (gravity is locally effaced in a freely 
falling reference frame) and the "strong" one (the gravitational mass of a system equals its inertial 
mass to which all forms of energy, including gravitational energy, contribute). 

The first principle holds because matter is minimally coupled to the metric of a curved spacetime 
so that test bodies follow geodesies. The second holds because minkowskian coordinates can be used 
in the vicinity of any event. The fact that the latter, strong, principle holds is ultimately due to 
the existence of superpotentials which allow to dehne the inertial mass of a gravitating system by 
means of its asymptotic gravitational field, that is, in terms of its gravitational mass. 

Nordstrom's theory of gravity, which describes gravity by a scalar field in flat spacetime, is 
observationally ruled out. ft is however the only theory of gravity with General Relativity to obey 
the strong equivalence principle. I show in this paper that this remarkable property is true beyond 
post-newtonian level and can be related to the existence of a "Nordstrom-Katz" superpotential. 

I. INTRODUCTION 

Among scalar-tensor theories of gravity, where gravity is described by a 4-dimensional metric g^ u together with a 
scalar field two stand apart, see e.g. [1] : General Relativity which describes gravity by a metric alone to which 
matter is minimally coupled, and Nordstrom's theory [2] (see [3] for an historical perspective) which describes it by a 
scalar field in flat spacetime ( "Einstein frame" formulation) . These two theories share the unique property to embody 
the strong equivalence principle, that is : the gravitational mass of a system is equal to its inertial mass, sum of the 
masses of its components and of all their binding, including gravitational, energies. 

The reason for which this principle holds in General Relativity is twofold : (1) the theory is a purely metric theory 
and there exist superpotentials, e.g. the Katz or KBL superpotential [4], which allow to define the inertial mass M- ln 
of a gravitating system by means of its gravitational field at infinity, that is, in terms of its (active) gravitational mass 
A/ g ; (2) the theory is second order. 1 

The claim that the principle holds in Nordstrom's theory as well (see e.g. [1]) relies on the fact that it can also 
be formulated as a purely, second order, metric theory ("Jordan frame" description). However this argument is not 
compelling since no superpotential for Nordstrom's theory has been proposed so far. 

I show here on various examples how, indeed, we do have M- ln = M g in Nordstrom's theory, even when the 
gravitational energy contributes significantly to M in (and not just at post-newtonian order). I shall also propose a 
superpotential for that theory, which, unfortunately, does not completely determine the conserved charge M- ln ... but 
the reason for this ambiguity seemed to me interesting enough to be given brief attention. 
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1 Indeed, metric theories derived for example from a lagrangian which is a function of the scalar curvature ("/(-R)" theories) or which 
are quadratic in the Riemann tensor do possess superpotentials, but do not obey the strong equivalence principle, see e.g. [5]. On the 
other hand, Gauss-Bonnet or Lovelock gravity theories, which are second-order and metric (but trivial in dimension four), do obey the 
strong equivalence principle, see e.g. [6]. 



II. NORDSTROM'S FIELD EQUATIONS IN THE EINSTEIN FRAME 



The action for gravity is taken to be 
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(2.1) 



c and G are the speed of light and Newton's constant ; the coefficients of the Minkowski metric and its determinant 
are £ M „ and I in the coordinate system x^ (and reduce to = r}^ = (— 1,+1,+1,+1) in an inertial frame with 
cartesian coordinates). The potential $(x M ) is dimensionless. 

If matter is an ensemble of particles with (inertial) mass m and proper velocities = = x^, the action S m 
describing its interaction with gravity and the corresponding stress-energy tensor T™ = — -j= are : 



i-^mc 2 J(l + F(<S>))dr , = ^mcj(l + F($))^^64x x -x x (T)])dT 
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where r is the proper time along their worldline, such that i^u^ u v = —c 2 and where is an a priori arbitrary 

function of $ (which must tend to 3> in the newtonian limit). If matter is a perfect fluid one takes : 



(1 + F($)) 4 (e+p) 
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(2.3) 



e and p being its energy density and pressure (the rationale for the various couplings to $ is given in [7] , see also [8] , 
and will become transparent in section 7 below). 

The equations of motion extremise (S g + S m ). They are : 



AttG dF/d® ^ 
c 4 1 + F 

1 + r 



(2.4) 



where D and □ are the covariant derivative and dalembertian associated with 



III. WEAK AND EINSTEIN'S EQUIVALENCE PRINCIPLES 

The fact that the coupling constant between a particle and gravity is its inertial mass m embodies the weak 
equivalence principle : all particles will fall the same way in a gravity field. In fact, the equations of motion for 
particles, (2.4b) with (2.2), can be rewritten as : 



~dV 
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where, as expected, m docs not appear. 

As an aside let us mention here that particles with zero mass travel at c along Minkowski's light cones and suffer no 
deviation ; hence the PPN parameter 7 is 7 = — 1 (instead of 7 = 1 as in General Relativity) . Let us also mention the 
value of the perihelion advance : expanding F($) as F($) = $ + ia2 ( E >2 + ... with $ = — — , a standard calculation 
yields Auj = — ^^Awgr where Acjqr is the general relativistic value, so that f3 — where (3 is the PPN parame- 
ter defined by 2 i~P+ 2 = ■ Nordstrom's theories are therefore observationally ruled out, whatever the function F. 

Showing that Nordstrom's theories obey Einstein's equivalence principle is a standard exercise. 
One first introduces Fermi coordinates: 

Let = cP(t), where r is proper time, be a worldline in some inertial frame with minkowskian coordinates 
X^. Suppose for simplicity that the motion is vertical : d M (r) = (d° (t) , d 3 (r j) . Go from the initial Lorentz frame 
(e ,e 3 ) attached to the origin O, to the tangent frame attached to a point O 1 on the worldline and defined as 

e' = d/c and e 3 = -d/VcU, (for d 3 < 0). Consider a point P. We have OP = OO' + O'P, that is, cTe + Ze 3 = 
d°e + d 3 e 3 + cT'e' + Z'e' 3 . For each P there exists, at least in a neighbourhood of the worldline, a unique point O' 
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and hence proper time r = t such that T" = 0. Setting Z' = z, the transformation from the minkowskian coordinates 
(T, Z) to the Fermi coordinates (t, z) is hence defined by 

Z = d 3 (t) + zy/l + (dP/c) 2 , cT = d a {t)+zd 3 /c. (3.2) 
In this coordinate system the equation of the worldine = cP(t) is z = and the metric becomes 

ds 2 = -c 2 dT 2 + dZ 2 = - ( l + z ^ 3/c2 | c 2 dt 2 +dz 2 . (3.3) 

V sJl + ^/c) 2 ) 

Impose now that the worldline = cP(t) is that of a test particle moving in a given gravitational field, that is 
that oP(t) solves (1) for a given function F($) and a given potential $(X M ). We then have that 



^i + ^/c) 2 



T=d".Z=d 3 ® Z 



(3.4) 

z=0 



and d° = y 1 + (d 3 /c) 2 . Hence the metric (3) becomes 

ds , = _ {1 _iit)z) 2 dt 2 + dz 2 where ^ 

az 



#• (3-5) 



Let us now consider another test particle with worldline t — t(r), z — z(t) where t is its proper time, whose 
equation of motion is also given by (1). In the Fermi coordinates (t,z), where now D is the covariant derivative 
associated with the metric (5), it reads : 



Z ~ l-zg/c 2 8t +(1 + Z /C '{l-zg/c 2 ° 8z ) ' *" 1 - gz/c 2 ' ^ 

When c — > oo these equations reduce to z = g(t) — c 2 ^ , t — t so that we recover the well-known newtonian result 
that the motion is uniform, z — and, thus, that gravity is effaced in the accelerated, Milne, frame, if the field is 
constrained to be uniform in z : = g(t)z/c 2 + $>o(t)- The equation of motion (4) for the origin O' of the frame is 
then given by d 3 = —g(t) in the original, inertial, frame. 

When zg/c 2 is no longer negligible, then all particles with zero initial velocities will remain at rest and gravity will 
be effaced in the Fermi frame if 1 _f g / c 2 = c 2 §7j that is, if ^ = — ln(l — g{t)z/c 2 ) + ^o{t) and, again, the equation 

of motion of the origin O' in the original, inertial frame, is given by (4) and reads -, d . = —g(t). If g — Const. 

y/l + (d 3 /c) 2 

the solution is uniformly accelerated motion. (For examples of motion in various gravitational fields, see e.g. [8].) 



IV. CONSERVATION LAW AND THE TOTAL MASS OF A GRAVITATING SYSTEM 

Let us now turn to the strong equivalence principle. 

The system being closed, it follows from the equations of motion (2.4) that the total stress-energy tensor is con- 
served : 

£>„(T£"+Tr)=0 with T^ = -^(d li *d v *-±t lu ,d p *&>*) (4.1) 

and T£ v given in (2.2) or (2.3). 

In an inertial frame and in cartesian coordinates, (1) becomes d v {T£ u + T£ v ) = 0. Integrating over all space we 
thus deduce that 

^L = ^^T§ i n i dS with M in c 2 = JdV{T™ + T™) (4.2) 

where t is the time in the inertial frame, where S is the 2-sphere at infinity , where n l is the unit vector pointing out 
of S, and where dS is the volume element of S : n l = (1, 0, 0) and dS = r 2 sin 9 d6 d<p in spherical coordinates. Since 
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the motion is confined, only (and not Tq?) contributes to the surface integral. In the definition of M- ln the first 
integral is over all space, the second is over the bodies creating the field. M- ln is, by definition, the total (inertial) 
mass of the system in the inertial frame where the 3-impulsion vanishes, J dV (T* + Tj^') = 0. 

Now, outside the source and far away, when the field is static or when radiation can be neglected, the field equation 
(2.4) for $ reduces A$ = with solution 

$ — > r-^ (4.3) 

c z r 

where r is the (large) distance from the source and M g the "active gravitational mass" of the system. In any specific 
problem, M g is related to either the mass to of the particles or the energy density e of the fluid creating the field, 
that is, ultimately, to M- m which, when the field is static or when radiation can be neglected, is then constant. 

The remarkable property of Nordstrom's theories is that, when F($) = $, then M- ln = M g and hence the theory 
obeys the strong equivalence principle. 2 

I will now show on specific examples how this result comes about. 



V. THE EXAMPLE OF THE TWO-BODY MOTION AT LOWEST ORDER 



Let us compute the gravitational field created by two point-like particles. 

In an inertial frame with cartesian coordinates = (t,r) the field equations are, see (2.2) and (2.4) : 

= ^ E m / dr % 6 ^ x " ■ ( 5J ) 

Expanding F as 

F($) = $ + ifl 2 $ 2 + - (5.2) 



this equation is solved iteratively, following the method set in [9] and [10], see also [11] and [13]. 

dF 



At lowest order, 4£ = 1 and the solution is the Lienard-Wiechert potential 



$(^) = -^|^+0(G 2 ) (5.3) 

d X* 1 

where, if is the intersection of the past light cone of with the worldline of to and w R = the 4-velocity of 
to at X£, then r R = -(X» - X£) UflR /c. 

At next order = 1 — ^ {r^ a 2 which, when evaluated on the wordline of to, is renormalized to jpg = 1 — , 

where, if X R is the intersection of the past light cone of X^ with the wordline of to' and u£ the 4-velocity of to' at 

X R , then /)r = — (X^ — X^)u' liR /c. Therefore the potential becomes 

= - E ^ + E + • ( 5 - 4 ) 

^ c 2 r R ^ c 4 r R p R 

The next steps are standard : perform a 3 + 1 decomposition (u R = c/^/l — v' R /c 2 etc) ; Taylor expand in 1/r 
where r = nr is the (large) separation between the point r and a reference point O in the system ; Taylor expand 
1/c. Express all quantities at time to = t — r/c. The final result is 



*(^)-e5?{i+^4 [MX*-*) + M) a - 1- 2 } } + e 



G 2 mm'a,2 
c 4 rR 



+ ■■■ (5.5) 



where z, v and v are the position, 3- velocity and acceleration of m at time to and where R is the distance between 
the two particles : RN = z '— z 1 '. 



2 There is however an exception, which was pointed out to me by Stanley Deser : an electromagnetic wave packet, whose stress-energy 
tensor is traceless, creates no gravitational field and hence has no gravitational mass; its inertial mass however is not zero since it is 
oc J dV (E 2 + B 2 ) in the frame where the Poynting vector is zero. 
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We then go to the center-of-mass inertial frame where the 3- impulsion of the system vanishes, z = m'RN/M + 
C(l/c 2 ) with M = m + m', so that the gravitational potential becomes (with V = R): 

. GM Gram!,, ■., . , T . 9 ~ G 2 mm'ai ,„ „. 

^) = -^M^ [{n - V){n - R) + {n - V) V /2] + 2 c ±rR + ' ' ' (5 ' 6) 

At the order considered here the motion is newtonian : V = —GMN/R and, for simplicity, we shall assume that 
the motion is circular so that V 2 = GM /R with R constant ; finally one takes the average on the orbital motion so 

that : (nV) 2 - ^{nNf = 0. 

Therefore, all in all, one obtains, at lowest order : 

GM 

$(r,t) = r-^ with M„ = M 

c 2 r s 



Gr 



2MRc 2 



(l + 4a 2 ) 



(5.7) 



Now, the inertial mass of the system is given by (4.2) with the stress-energy tensors given by (4.1) and (2.2). At 
the order considered one has 

J d VT™ = Y, mc2 ( 1 + ^+*)+--- > |^T g 00 = -^m$ + ... (5.8) 

where $ is given by (3), that is, at lowest order, by $ = J2 §7 and must be evaluated on the (newtonian) trajectories. 
After renormalization we have, in the center-of-mass frame and for a circular orbit : 

, , 9 \ -« 9 1 \ ^ 9 Gram! . , 9 /, Gram! \ . . 

M in c 2 =Y.^ + - 2 Y.™* - — = Mc 2 (l - — 2 j (5.9) 

and hence 

^ = 1 + ^ with V = 2a 2 (5.10) 

where E g — — Gm ^ a is the (Newtonian) gravitational energy of the system and 77 Nordtvedt's PPN parameter. (One 
notes, that as in all scalar-tensor theories, see [1], r\ is related to (3 and 7 by the PPN relation : r\ = 4/3 — 7 — 3.) 

We therefore see on this example that the (active) gravitational mass M g which appears in the potential far away 
from the system, will be equal to the inertial mass M- m of the system if a-i — 0. 3 

One could of course proceed and thus constrain the function F(Q), order by order, but the calculations soon become 
heavy (although much simpler than in General Relativity, see eg [11]) and will not show, in any case, what the function 
F should be for the strong equivalence principle to hold exactly. This is why we proceed to the next example. 4 



VI. THE EXAMPLE OF A CONSTANT DENSITY "STAR" 



If the gravitational field is static (thus guaranteeing the absence of radiation), the integral (4.2) giving M- m can be 
reduced to an integral over the body creating the field using the field equation (2.4) and taking into account that 
$ oc i at infinity : 



3 If the system is electromagnetically, rather than gravitationally bound, then the gravitational field it creates is also given by (1- 
7) since the electromagnetic stress-energy tensor, being traceless, does not contribute. At lowest order the motion is coulombian, 

V = qq' MN /mm' R 2 so that, for 02 = : M g = M + qq' /2Rc 2 , which, again, is its inertial mass. 

4 As an aside and a small tribute to Joshua's seminal contribution to the problem of motion in General Relativity [9], I give in the 
Appendix the gravitational energy lost by a binary system in Nordstrom's theories. 
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Therefore, for matter being a perfect fluid with stress-energy tensor given in (2.3) : 



M in c A = 



dV(l + F) 3 



where the equation fo $ is 



47rG._ ^*dF . 
A$ = - — ( 1 + F ) 3 ^(^- £ ) 



and the Euler equation for matter reduces to 



e + p 



dF/d<5> 
' l + F 



(6.2) 



(6.3) 



(6.4) 



If, moreover, the configuration is spherically symmetric so that $, e and p depend on the radial coordinate r only, 
then A<£> = ^-(r 2 (f>')', Vp — > p', with a prime denoting derivation with respect to r, and the solution for $ outside the 

The junction conditions are that p be zero at the surface of the star and $ and $' be continuous. 



star is $ 



An equation of state is required to close the system of equations. We shall consider here as an example the unrealistic 
but simple case when e = Const (the more realistic case of a barotropic fluid is studied numerically in [12]). Then 
the Euler equation (4) with the condition that p vanishes at the surface of the star integrates as 



P - 



F R -F 
l + F 



where F R = R being the radius of the star. As for the equation (3) for $ it becomes 



^*y = ^(l + Ff%(l-3F R + 4F). 



(6.5) 



(6.6) 



Outside the star : $ = with $ fl 

its derivative are 



- G ]f R an integration constant, so that the conditions of continuity of <£ and 

(6.7) 



+ m\R) = , $(R)-$ R = 0. 
Once the solution of (6-7) for $(r) is obtained the inertial mass (2) is given by 



M in c 2 = 47re / drr 2 (l + F) 2 
Jo 



(1 + F)2 ~2 $ d¥ (1 ~ 3F * + 4 ^ 



(6.8) 



Let us start with the post-newtonian approximation. Expanding F($) as F(<&) = $ + ^a 2 < E> 2 + • • • it is an exercise 
to solve (6) iteratively with the boundary conditions (7) to obtain 



$ = 



3$ fl AwGe r 2 



and 



AirGe 



R 2 



-3$; 



1_ ?*« ( 7 + 20,) 



(6.9) 



with <P R = -^f . At the order considered (8) reduces to M in c 2 = Aire f*(l + ^) so that M in = M g (l - 6*^2) 
and we recover on this example the result obtained in the previous section, that is 



— - = 1 + n — — 

M- m 'Mc 2 



with n = 2a 2 



(6.10) 



where E g — — 3( ^^ is the newtonian gravitational energy of the body and where 77 is the Nordtvedt parameter. 



5 A word of caution : from the conservation law (4.1) one can deduce that 5^3- JdVx'x^T 00 = JdVT'i up to surface terms. In the 
static and spherically symmetric case considered here one may therefore be tempted to conclude that J dV T lJ = 0, so that Jdr r 2 T£f = 
- fdrr 2 T^ r . Now since 7J r = T°° this would imply that M in c 2 = AnJdrr 2 (T^> - T£~), that is : M in c 2 = 4tt fdrr 2 (l + F) 2 (e-p) 
when matter is a perfect fluid. For constant e, this formula, using (5) and (9) gives a wrong result, different from (10), the reason being 
that the integral fdVx''x : >T 00 does not converge. 
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Now, it is easy to go beyond the post-newtonian approximation as equation (6) can be solved numerically : for any 
specific function F($) and value for one chooses an initial value for $ and to wit : $(0) = $o and $'(0) = 0, 
and fits $o so that after integration up to some R, the junction conditions (7) are satisfied at R. 

Once the solution for <f>(r) is thus obtained one can compute the value of the inertial mass (8). If M- ln = M g or 
equivalently, —GM g /c 2 R — then we must find that, for all 3>h : 



R*R + ^ l R drr\l + Ff 
c Jo 



1 HF 

(1 + F)2 ~2 $ rf¥ (1 ~ 3 ^ + 4F) 



0. (6.11) 



The result (using e.g. Mathematica) is that this is indeed true if the function F($) is 

F($) = $. (6.12) 

We have therefore shown that the strong equivalence principle holds exactly (in this particular case of a constant 
density "star") in Nordstrom's "final" theory [1] [2] where F($) = One could look at other examples with more 
realistic equations of state (e.g. polytropic as in [12]) but, again, this would not prove that the strong equivalence 
principle holds exactly in all cases. This is why we turn now to the Jordan frame description of Nordstrom's theories. 



VII. NORDSTROM'S FIELD EQUATIONS IN THE JORDAN FRAME 

As was already known to Einstein [2-3] Nordstrom's final theory of gravity can be turned into a metric theory. 
Indeed if one sets 

g lu , = (l + F) 2 e iu , (7.1) 
then the equations of motion (2.4) can be recast as 

24ttG fdF\ 2 - 6 d 2 F, 



R =—[m) T --TTF^W^ (7.2) 

I A/fr = o 

where R = — ^ 1+ 6 F ^ DF is the scalar curvature of the conformally fiat metric g pv and where T M „ is the stress-energy 
tensor of matter minimally coupled to the metric g^ v . Thus the action and stress-energy tensor for particles given in 
(2.2) become 

S m = -^mc 2 Jdf , TT = ^mc J ^^S 4 (x x - x x {f))df (7.3) 
with {t M = and g^u^u" = — c? . As for the stress-energy tensor for a perfect fluid (2.4) it reads : 

f™=(e+p)^+p~g,„. (7.4) 

In this "Jordan frame" formulation the special status of Nordstrom's final theory, F($) = <&, jumps to the eye. In 
that case indeed the equations of the theory reduce to 

R=^-T m , D v f^=0 , C^ pa = (7.5) 

where the vanishing of the Weyl tensor C^ pa imposes the metric to be conformally flat. Equations (5) share with 
Einstein's equations the fact that they are purely geometrical and second order. Hence the claim, cf e.g. [1], that 
Nordstrom's final theory embodies the strong equivalence principle. 

The problem however is that the conservation law (4.1) translates into D V T^ = and, since spacetime is no longer 
flat, there is no coordinate system which reduces it to d^T^ = 0. Therefore, just as in General Relativity, there is 
no obvious conservation law from which to compute the inertial mass of a system. In order to find one the action, 
functional of the Jordan frame metric g^ must be found which gives the equations of motion (5). 
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VIII. A JORDAN ACTION FOR NORDSTROM'S GRAVITY 



N.B. : I shall henceforth drop the tildes which decorate the formulas of the previous section. 
Let us consider the following action : 6 

where g is the determinant of the (Jordan) metric g^ v , where R and C^ vpa are the corresponding scalar curvature 
and Weyl tensor and where \p Vpa is a Lagrange multiplier possessing all the symmetries of the Weyl tensor. Note 
that the first term is minus one third the Einstein-Hilbert action for General Relativity. As for matter we take it to 
be minimally coupled to the metric g p „ so that its action S m is that of Special Relativity with — > g^ v . 
Extremization with respect to the matter variables gives, as usual 

D V T» V = (8.2) 

where T^ v = — -j=jMj^ is its stress-energy tensor. Extremization with respect to \^ pa imposes 

C^vpa = . (8-3) 

Finally, extremization with respect to g^ v (ignoring boundary terms for the time being) yields an equation whose 
trace is 

R=^T m (8.4) 



and whose traceless part is 



2D a D^J + K$ \J = { - \9„tA (r, v \g, v R ) . ( K r, , 



Equations (2-4) are the Jordan frame version of Nordstrom's final theory, see (7.5) : Eq. (3) imposes the metric to 
be conformally flat, g pv = (1 + ( f>) 2 ^ AlJ , ; Eq. (2) and (4) can then be recast in their original, Einstein frame, version 
Eq (2.4) (with F = $) and determine the conformal factor $(x M ) and the motion of matter. 

As for (5) it can be rewritten in terms of the flat metric and its covariant derivative that we now ornate with 
a bar D. Indeed we have that 

2D a D p \ IM J + R1\J = ^f^/ , (8-6) 

i? ^~" 2 TT^"^TT^ + 4 (iTW _V (rTW ' JTTW (8J) 

Therefore, once the solution for $ is known, e.g., $ = —-^f- outside a static and spherically symmetric distribution, 
then (5) together with (6-7) is an equation for the Lagrange multiplier \y, vpa . However, \^ vpa , having the symmetries 
of the Weyl tensor, possesses ten independent components, whereas (5), being traceless, has only nine components. 
The system of equations (5) for A M „ p(T is therefore undetermined. This is of no consequence to obtain the gravitational 
field <& since \p Vpa docs not enter its equations of motion. But this under-determination will prevent the action (1) 
to yield a well defined inertial mass as we shall now see. 

IX. A JORDAN FRAME DEFINITION OF INERTIAL MASS IN NORDSTROM'S THEORY 

A. Katz superpotential and conserved charges 

I give here a brief (and hopefully comprehensible) summary of how to build a superpotential out of a metric 
lagrangian for gravity. For details see [4] and [6] . 



I thank Misao Sasaki for suggesting it. 
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Consider the lagrangian density ^fj, k being some coupling constant, with 

C=y/^gC and C = L + D^ (9.1) 

where L is a scalar, functional of the metric g pv and its derivatives up to the second, and where k p is some vector. 
Its variation with respect to the metric can be written as 

Si = -a^5g„ v + d^W + 6hP) (9.2) 

with 

V^ = a^"5g vp + ^P a 5TZ p (9.3) 

where Y a vp are the Christoffcl symbols and where a p „, a A " yp and fi tlvp a are some tensors depending on the specific form 
of the lagrangian L. 

If, now, the variation Hg pv is due to a mere change of coordinates then 6 reduces to a Lie derivative and it is an 
exercise to see that (2) can be cast into the following form : 

8^ = 2^0^ (9.4) 
where the "current" j p = j£ + j£ is given by (parentheses denoting symmetrization, brackets antisymmetrization) 

f a = (i<r + + R\ pa ) - 2a^ D v i p - (3^^D^ a 

Now the right hand side of (4) is identically zero by virtue of the (generalized) Bianchi identity. Therefore the current 
is identically conserved : d^j^ = . 

The conservation of j p implies that there exists an antisymmetric "superpotential" such that 

f = D v jW = dJjW . (9.6) 
Looking for an expression of the form = + with 

j.JH = A l^]p^ p + B^ pa D p ^ and j^ v] = 2^k^ (9.7) 

we get from (5) 

Having thus constructed j^ v \ the Katz superpotential and charge are defined as 

Q = --^ [ d 2 x nJW with = jW _ . (9.9) 

2kc z J s 

In this formula, S is the 2-sphere at infinity and d 2 x — sm9d9d(j) with = (1,0,0) in asymptotically spherical 

coordinates. As for it is the superpotential corresponding to another (background) metric g pv and serves as a 
regulator. 

The charge Q depends on the vector k p , which is chosen in order that the boundary conditions of the variational 
principle be Dirichlet's. It also depends on the vector : we shall choose it to be the Killing vector corresponding 
to time translations at infinity so that Q is then the inertial mass of the system. 

B. Application to Nordstrom's theory 

Let us apply now this machinery when the Lagrangian L in (1) is that of Nordstrom's gravity, that is, 

L = R + \;r°C» vpa . (9.10) 
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Computing its variational derivative with respect to the metric yields a vector V 1 as given in (2) with 

a ti{vp) _ — 2Z? cr A M ^ p ^ <T , pv( v p)< 7 — g' /p g fta — g^ v g p ^ a — 2\^ l/p "' a . (9-H) 
Solving (8) then gives : 

ft\MP° = gUPg"" _ gVPgH* + , ^IHP = AD ' a ' P ° . (9.12) 

Therefore the superpotential and charge are defined in (9) with, cf (7) and (12) : 

f v = 2D^"] + 2\^ pa D p £ a + ^ p D a \ pvpa + 2£^ . (9.13) 

As for the vector fc p we choose it so that its variation <5fc p cancels out the terms proportional to 5Y° p in the divergence 
(2.3): 

k» = -fl pvp a Af p = fA; p - g p °A% + 2A^A^ where = 1^ - f ^ . (9.14) 

In the case of Einstein's theory where the \^ pa terms are absent, j pu and fc M as given in (13) (14) are the 
superpotential and vector first proposed in [4]. 7 



C. Nordstrom-Katz' inertial mass 



In Nordstrom's theory the (Jordan) metric outside a spherically symmetric body is 

ds 2 = {I + <t>) 2 [-dt 2 + dr 2 +r 2 (d0 2 + sin 2 6 d<j> 2 )] with * = -^ i - (9-15) 

Since this metric is not Schwarzschild's and since the coupling constant is minus one third that of Einstein's, 
« = _24^g S ee (8.1), the contribution of the Einstein part of the lagrangian (10) to the inertial mass (that is for 
£ p = (1,0,0,0)) turns out to be plus one third of the gravitational mass M g . 

Using now the symmetries of A A " ypCT and the expression of the Christoffel symbols for the metric (15), 

r -p = f -p + TT¥ ( ^ 9p$ + s p d ^ - LpB ^ (9A6) 

where f p p are the Christoffel symbols of the flat background metric £ M „ in spherical coordinates, it is easy to see that 
the term 2A M!/p J A^ in the expression (14) of the vector fc M as well as the term 2\ pvpa D p ^ a in the expression (13) of 
do not contribute. Therefore, all in all, we have 

3 3G r- 

We thus see that, if the Lagrange field \ pupa plays no role in determining the gravitational field, it docs enter the 
definition of the inertial mass of a gravitating system. 

Now, as we have seen in section 8, the equations of motion (8.5-7) for A A " ypCT do not determine it completely. A 
closer look at these equations when the metric is given by (15) tells us that 

lim r 2 D a \ 0r0 ° (9.18) 
cr 

M in = ^(l-a). (9.19) 

The only way that I see to fix the value of the constant a is to return to the Einstein frame formulation of 
Nordstrom's theory and use the result obtained in section 6, to wit that, indeed, M; n = M g in the particular case of 
a constant density perfect fluid star. In this case then a — — 2 and, since the result (19) is general, then a must be 
equal to —2 in all cases. 



M in = - — Jim r 2 D a X 0r0 a . (9.17) 



with a an arbitrary constant. Therefore 



7 The metric outside a spherically symmetric object being Schwarzschild's in GR, ds 2 = — (l — 222) dt 2 + 1 _ d 2 m /r r ' 2 {dQ 2 + sm2 6 dtji 2 ) 

with m = G ^i S where M g is its gravitational mass, the computation of the inertial mass of the body, with £ M = (1,0,0,0), k = 
and a fiat background, gives Q = M[ n = M g . 

Now, the Schwarzschild solution also solves the field equations of pure quadratic theories L = aR 2 + fiR lJlv W J - v . However, as a 
straightforward generalization of the results presented in [6] shows, the inertial mass, as defined in (9), then vanishes, in agreement with 
[5] : the strong equivalence principle is violated in such theories. 
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X. CONCLUSION 



In this paper I explored some aspects of Nordstrom's theories of gravity which are of some interest, not only in an 
historical perspective, see [3] , but also because they shed some light on the thorny issue of the validity of the strong 
equivalence principle in relativistic theories of gravity. 

I first showed on the simple example of a perfect fluid constant density star, that, when formulated within a Special 
Relativity framework, Nordstrom's "final" theory does satisfy the strong equivalence principle exactly and not only at 
post-newtonian order. I then tried to show that that was always true by giving a metric, Jordan frame, formulation 
of the theory. I partly succeeded by exhibiting a Katz superpotential and associated inertial mass of a gravitating 
system which is indeed proportional to its gravitational mass. However I find that my argumentation to claim that, 
hence, the strong equivalence principle is always true in Nordstrom's final theory is a bit weak... 

A way to straighten the proof would be to start from a Jordan action different from the one I introduced in section 
8. A candidate, inspired by Ref. [14], could be, instead of (8.1), the Palatini-like action : 

S s [^,T^,\ a ^,f^} = -^ j d^il + ^f^gg^R^ - A- J d ^ X ^gX^ S R^ 5 (10.1) 

with 

= dpT^ — d v Y p Pil + r^r^ — r^r^p , (10.2) 

and 

R^s = d-y^ps ~ % r ^ 7 + r 7P r /35 ~ r ?p r ia 7 ■ (10.3) 

As for the action for matter it would be its Einstein frame version, see e.g. (2.2). 

As can easily be seen, extremization with respect to the (independent) connexion T^ 7 , with respect to $ and with 

respect to \£ jS yields the equations of motion (8.2-4) of Nordstrom's final theory. Finally extremization with respect 
to to f a &1 gives : 

D a \» aa = . (10.4) 

The next step would be to build a "Nordstrom-Katz" superpotential out of the above Palatini action. In order to 
do so the techniques developped in [15] will have to be used. This is left to further work. 
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Appendix A: Energy loss of a binary system 



When the system is radiating the conservation law (4.1) together with the definition (4.2) implies that the system 
looses energy at a rate given by 



dM h 
dt 



nyds 



with 



rpg. _ '- 

oi ~ AttG 



(Al) 



where t is the time in some inertial frame, where S is the 2-sphere at infinity , where n % is the unit vector pointing 
out of S and where dS is the volume element of S : n l = (1,0,0) and dS = r 2 sm9d9d<p in spherical coordinates. 
Since the motion is confined, only (and not T^) contributes to the surface integral. 

The gravitational potential created by two point-like bodies is given by (5.4) at C(G 3 ) but the Taylor expansion in 
1/c must now be pushed one order beyond that obtained in (5.5). The result is 



$(r,*) = $ [2] (f, 



1 \ ^ G mm'ao 



c 4 ri? 



(n.v) 



+ 3(n.v)(n.z){n.v) 
(N.V)(n.z)~ 



(n.z)(v.v) + (n.v) 3 — v 2 (n.v) 



R 



(A2) 
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where $r 2 ](r, t) is the lowest order expansion of $ given in (5.5). 

The next steps are : compute the derivatives do<& and <9;$ in the (newtonian) center of mass frame where z 
m'RN/M ; compute using the relations J riiUj 



dQ = etc ; use the newtonian equations of motion, V = 



-GMN I R to obtain, at the end of the day : 



dM h 
dt 



G fd 3 QY 
{dt*) 



4G 3 (t 



,^2 



30c 



9 c 5 R 4 



-(N.V) 2 (l + 3a 2 ) 2 



(A3) 



where Qij = ^^j-R 2 (3NiNj — Sij) is the quadrupole moment of the system. The first term is 1/6 the general 
relativistic value. Finally, since the orbit is keplerian at this order (RN — [a(cos?7 — e) , ayl — e 2 sin tj], t = 
\J gTT ~ e sin ^7) ) where a and e are its semi-major axis and eccentricity), one can compute the average loss over one 
orbit, AMjn = f Q P ^f^dt where P = 2irJ ^jg, and rewrite (3) as 



P 



1 AM, 
6" 



P 

AMj, 
P 



I GR 



I GR 



5 e 2 (l + ie 2 )(l + 3a 2 ) 5 



24 (1 + 1^ + -^ 



73 „2 



32 G 4 (mm') 2 M (l + ge 



96 c J 



(A4) 



(1 - e 2 ) 7 /2 



These results extend those obtained in [1] and [12] and may serve as a benchmark for testing numerical codes, see e.g. 
[12] [13]. 
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